3. Quantum Gases 



In this section we will discuss situations where quantum effects are important. We'll still 
restrict attention to gases — meaning a bunch of particles moving around and barely 
interacting — but one of the first things we'll see is how versatile the idea of a gas 
can be in the quantum world. We'll use it to understand not just the traditional gases 
that we met in the previous section but also light and, ironically, certain properties of 
solids. In the latter part of this section, we will look at what happens to gases at low 
temperatures where their behaviour is dominated by quantum statistics. 

3.1 Density of States 

We start by introducing the important concept of the density of states. To illustrate 
this, we'll return once again to the ideal gas trapped in a box with sides of length 
L and volume V = L 3 . Viewed quantum mechanically, each particle is described 
by a wavefunction. We'll impose periodic boundary conditions on this wavefunction 
(although none of the physics that we'll discuss in this course will be sensitive to the 
choice of boundary condition). If there are no interactions between particles, the energy 
eigenstates are simply plane waves, 
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Boundary conditions require that the wavevector k = (hi, k 2 , k 3 ) is quantized as 
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and the energy of the particle is 
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with k — The quantum mechanical single particle partition function (1.21) is given 
by the sum over all energy eigenstates, 



The question is: how do we do the sum? The simplest way is to approximate it by an 
integral. Recall from the previous section that the thermal wavelength of the particle 
is defined to be 
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The exponents that appear in the sum are all of the form ~ X 2 n 2 /L 2 , up to some 
constant factors. For any macroscopic size box, A C L (a serious understatement! 
Actually A <^gC L) which ensures that there are many states with < k B T all of which 
contribute to the sum. (There will be an exception to this at very low temperatures 
which will be the focus of Section 3.5.3). We therefore lose very little by approximating 
the sum by an integral. We can write the measure of this integral as 

£ x I d3n = wrl d3k= w?[ dkk? 

where, in the last equality, we have integrated over the angular directions to get 4tt, the 
area of the 2-sphere, leaving an integration over the magnitude k = \k\ and the Jacobian 
factor k 2 . For future applications, it will prove more useful to change integration 
variables at this stage. We work instead with the energy, 
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We can now write out integral as 
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where 

is the density of states: g(E)dE counts the number of states with energy between E 
and E + dE. Notice that we haven't actually done the integral over E in (3.1); instead 
this is to be viewed as a measure which we can integrate over any function f(E) of our 
choosing. 

There is nothing particularly quantum mechanical about the density of states. In- 
deed, in the derivation above we have replaced the quantum sum with an integral over 
momenta which actually looks rather classical. Nonetheless, as we encounter more and 
more different types of gases, we'll see that the density of states appears in all the 
calculations and it is a useful quantity to have at our disposal. 

3.1.1 Relativistic Systems 

Relativistic particles moving in d — 3 + 1 spacetime dimensions have kinetic energy 

E = Vh 2 k 2 c 2 + m 2 c 4 (3.3) 
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Repeating the steps above, we find the density of states is given by 

9{E) = 2^ VE2 - ^ (3 ' 4) 
In particular, for massless particles, the density of states is 

m = ^ (3-5) 
3.2 Photons: Blackbody Radiation 

"It was an act of desperation. For six years I had struggled with the black- 
body theory. I knew the problem was fundamental and I knew the answer. I 
had to find a theoretical explanation at any cost, except for the inviolability 
of the two laws of thermodynamics" 

Max Planck 

We now turn to our first truly quantum gas: light. We will consider a gas of 
photons — the quanta of the electromagnetic field — and determine a number of its 
properties, including the distribution of wavelengths. Or, in other words, its colour. 

Below we will describe the colour of light at a fixed temperature. But this also applies 
(with a caveat) to the colour of any object at the same temperature. The argument for 
this is as follows: consider bathing the object inside the gas of photons. In equilibrium, 
the object sits at the same temperature as the photons, emitting as many photons as 
it absorbs. The colour of the object will therefore mimic that of the surrounding light. 

For a topic that's all about colour, a gas of photons is usually given a rather bland 
name — blackbody radiation. The reason for this is that any real object will exhibit 
absorption and emission lines due to its particular atomic make-up (this is the caveat 
mentioned above). We're not interested in these details; we only wish to compute the 
spectrum of photons that a body emits because it's hot. For this reason, one sometimes 
talks about an idealised body that absorbs photons of any wavelength and reflects none. 
At zero temperature, such an object would appear black: this is the blackbody of the 
title. We would like to understand its colour as we turn up the heat. 

To begin, we need some facts about photons. The energy of a photon is determined 
by its wavelength A or, equivalently, by its frequency u> = 2nc/\ to be 

E = huo 

This is a special case of the relativistic energy formula (3.3) for massless particles, 
m = 0. The frequency is related to the (magnitude of the) wavevector by u = kc. 



-63- 



Photons have two polarization states (one for each dimension transverse to the di- 
rection of propagation). To account for this, the density of states (3.5) should be 
multiplied by a factor of two. The number of states available to a single photon with 
energy between E and E + dE is therefore 

Equivalently, the number of states available to a single photon with frequency between 
ui and ui + du is 



where we've indulged in a slight abuse of notation since g{uS) is not the same function 
as g(E) but is instead defined by the equation above. It is also worth pointing out an 
easy mistake to make when performing these kinds of manipulations with densities of 
states: you need to remember to rescale the interval dE to du. This is most simply 
achieved by writing g(E)dE = g(uj)du as we have above. If you miss this then you'll 
get g(u>) wrong by a factor of H. 

The final fact that we need is important: photons are not conserved. If you put 
six atoms in a box then they will still be there when you come back a month later. 
This isn't true for photons. There's no reason that the walls of the box can't absorb 
one photon and then emit two. The number of photons in the world is not fixed. To 
demonstrate this, you simply need to turn off the light. 

Because photon number is not conserved, we're unable to define a chemical potential 
for photons. Indeed, even in the canonical ensemble we must already sum over states 
with different numbers of photons because these are all "accessible states". (It is 
sometimes stated that we should work in the grand canonical ensemble at \i = which 
is basically the same thing). This means that we should consider states with any 
number N of photons. 

We'll start by looking at photons with a definite frequency uj. A state with N such 
photons has energy E = Nftw. Summing over all N gives us the partition function for 
photons at fixed frequency, 



g(E)dE = g(uj)duj 
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We now need to sum over all possible frequencies. As we've seen a number of times, 
independent partition functions multiply, which means that the logs add. We only need 




Figure 14: The Planck Distribution function (Source: E. Schubert, Light Emitting Diodes). 

to know how many photon states there are with some frequency u. But this is what 
the density of states (3.6) tells us. We have 

/*oo t 7 roc 

\ogZ= / duj g{w) \ogZ w = — / duj oj 2 log (l - e - ^) (3.8) 
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3.2.1 Planck Distribution 

From the partition function (3.8) we can calculate all interesting quantities for a gas of 
light. For example, the energy density stored in the photon gas is 

E --lH^ osZ = ^L d "^T (3 ' 9) 

However, before we do the integral over frequency, there's some important information 
contained in the integrand itself: it tells us the amount of energy carried by photons 
with frequency between u and oj + du 

E (") d «= ( 3 - 10 ) 

This is the Planck distribution. It is plotted above for various temperatures. As you 
can see from the graph, for hot gases the maximum in the distribution occurs at a lower 
wavelength or, equivalently, at a higher frequency. We can easily determine where this 
maximum occurs by finding the solution to dE(u)/du = 0. It is 

n 

where ( ~ 2.822 solves 3 — ( = 3e~~^. The equation above is often called Wien's 
displacement law. Roughly speaking, it tells you the colour of a hot object. 
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To compute the total energy in the gas of photons, we need to do the integration 
in (3.9). To highlight how the total energy depends on temperature, it is useful to 
perform the rescaling x = flhu, to get 

V {k B T) 4 f°° x 3 dx 
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The integral I = J dxx 3 /(e x — 1) is tricky but doable. It turns out to be / = 7r 4 /15. 
(We will effectively prove this fact later in the course when we consider a more general 
class of integrals (3.27) which can be manipulated into the sum (3.28). The net result of 
this is to express the integral / above in terms of the Gamma function and the Riemann 
zeta function: I = r(4)£(4) = 7r 4 /15). We learn that the energy density £ = E/V in a 
gas of photons scales is proportional to T 4 , 

7T 2 t 4 
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Stefan-Boltzmann Law 

The expression for the energy density above is closely related to the Stefan-Boltzmann 
law which describes the energy emitted by an object at temperature T. That energy 
flux is defined as the rate of transfer of energy from the surface per unit area. It is 
given by 

Sc 

Energy Flux = — = aT A (3.11) 

where 

- - ^ k %. = 5.67 x 10^8 Js^m' 2 ^ 



Q0ft 3 c 2 
is the Stefan constant. 

The factor of the speed of light in the middle equation of (3.11) appears because 
the flux is the rate of transfer of energy. The factor of 1/4 comes because we're not 
considering the flux emitted by a point source, but rather by an actual object whose 
size is bigger than the wavelength of individual photons. This means that the photon 
are only emitted in one direction: away from the object, not into it. Moreover, we only 
care about the velocity perpendicular to the object, which is (ccos^) where 6 is the 
angle the photon makes with the normal. This means that rather than filling out a 
sphere of area 4n surrounding the object, the actual flux of photons from any point on 
the object's surface is given by 
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Radiation Pressure and Other Stuff 

All other quantities of interest can be computed from the free energy, 

F = -k R T log Z 
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We can remove the logarithm through an integration by parts to get, 
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From this we can compute the pressure due to electromagnetic radiation, 



dF 

P= ~dV 



W 3c 



This is the equation of state for a gas of photons. The middle equation tells us that the 
pressure of photons is one third of the energy density — a fact which will be important 
in the Cosmology course. 

We can also calculate the entropy S and heat capacity Cy. They are both most 
conveniently expressed in terms of the Stefan constant which hides most of the annoying 
factors, 
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3.2.2 The Cosmic Microwave Background Radiation 

The cosmic microwave background, or CMB, is the afterglow of the big bang, a uni- 
form light that fills the Universe. The intensity of this light was measured accurately by 
the FIRAS (far infrared absolute spectrophotometer) instrument on the COBE satellite 
in the early 1990s. The result is shown on the right, together with the theoretical curve 
for a blackbody spectrum at T = 2.725 K. It may look as if the error bars are large, 
but this is only because they have been multiplied by a factor of 400. If the error bars 
were drawn at the correct size, you wouldn't be able to to see them. 
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This result is totally astonishing. The light 
has been traveling for 13.7 billion years, almost 
since the beginning of time itself. And yet we can 
understand it with ridiculous accuracy using such 
a simple calculation. If you're not awed by this 
graph then you have no soul. 

If you peer closely enough, the data does depart 
from the blackbody curve. You need to measure 
to an accuracy of one part in 10 5 . That's where 
the story gets really interesting. But it's a topic 
for another course. 

3.2.3 The Birth of Quantum Mechanics 

The key conceptual input in the derivation of Planck's formula (3.10) is that light of 
frequency u comes in packets of energy E = Hlj. Historically, this was the first time 
that the idea of quanta arose in theoretical physics. 

Let's see what would happen in a classical world where light of frequency u can have 
arbitrarily low intensity and, correspondingly, arbitrarily low energy. This is effectively 
what happens in the regime hu <C k^T of the Planck distribution where the minimum 
energy hoj is completely swamped by the temperature. There we can approximate 

1 _ 1 

- 1 ~ ]3frw 

and Planck's distribution formula (3.10) reduces to 

Notice that all hints of the quantum h have vanished. This is the Rayleigh- Jeans law for 
the distribution of classical radiation. It has a serious problem if we try to extrapolate it 
to high frequencies since the total energy, E = J °° E(u)du, diverges. This was referred 
to as the ultra-violet catastrophe. In contrast, in Planck's formula (3.10) there is an 
exponential suppression at high frequencies. This arises because when hoo ^> ksT, the 
temperature is not high enough to create even a single photon. By imposing a minimum 
energy on these high frequency modes, quantum mechanics ensures that they become 
frozen out. 
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Figure 15: 



- 68 - 



3.2.4 Max Planck (1858-1947) 

"A new scientific truth does not triumph by convincing its opponents and 
making them see the light, but rather because its opponents eventually die" 

Max Planck 

Planck was educated in Munich and, after a brief spell in Kiel, moved to take a 
professorship in Berlin. (The same position that Boltzmann had failed to turn up for). 

For much of his early career, Planck was adamantly against the idea of atoms. In 
his view the laws of thermodynamics were inviolable. He could not accept that they 
sat on a foundation of probability. In 1882, he wrote "atomic theory, despite its great 
successes, will ultimately have to be abandoned" . 

Twenty years later, Planck had changed his mind. In 1900, he applied Boltzmann's 
statistical methods to photons to arrive at the result we derived above, providing the 
first hints of quantum mechanics. However, the key element of the derivation — that 
light comes in quanta — was not emphasised by Planck. Later, when quantum theory 
was developed, Planck refused once more to accept the idea of probability underlying 
physics. This time he did not change his mind. 

3.3 Phonons 

It is hard to imagine substances much more different than a gas 
and a solid. It is therefore quite surprising that we can employ 
our study of gases to accurately understand certain properties 
of solids. Consider a crystal of atoms of the type shown in the 
figure. The individual atoms are stuck fast in position: they 
certainly don't act like a gas. But the vibrations of the atoms 
— in other words, sound waves — can be treated using the 
same formalism that we introduced for photons. 

3.3.1 The Debye Model 

Quantum mechanics turns electromagnetic waves into discrete packets of energy called 
photons. In exactly the same way, sounds waves in solids also come in discrete packets. 
They are called phonons. We'll take the energy of a phonon to again be of the form 

E = Hlu = hkc s (3.12) 




where c s is now the speed of sound rather than the speed of light. 



The density of states for phonons is the same as that of photons (3.6) with two 
exceptions: we must replace the speed of light c with the speed of sound c s ; and phonons 
have three polarization states rather than two. There are two transverse polarizations 
(like the photon) but also a longitudinal mode. The density of states is therefore 



There is one further important difference between phonons and photons. While light 
waves can have arbitrarily high frequency, sound waves cannot. This is because high 
frequency waves have small wavelengths, A = 2ttc s /uj. But there is a minimum wave- 
length in the problem set by the spacing between atoms, it is not possible for sound 
waves to propagate through a solid with wavelength smaller than the atomic spacing 
because there's nothing in the middle there to shake. 

We will denote the maximum allowed phonon frequency aswjj. The minimum wave- 
length, \ D , should be somewhere around the lattice spacing between atoms, which is 
(V/N) 1 / 3 , so we expect that ujd ~ {N/VY^c s . But how can we work out the coeffi- 
cient? There is a rather clever argument to determine uod due to Debye. (So clever that 
he gets his initial on the frequency and his name on the section heading). We start by 
counting the number of single phonon states, 



The clever bit of the argument is to identify this with the number of degrees of freedom 
in the solid. This isn't immediately obvious. The number of degrees of freedom in a 
lattice of N atoms is 3N since each atom can move in three directions. But the number 
of single phonon states is counting collective vibrations of all the atoms together. Why 
should they be equal? 

To really see what's going on, one should compute the correct energy eigenstates 
of the lattice and just count the number of single phonon modes with wavevectors 
inside the first Brillouin zone. (You will learn about Brillouin zones in the Applications 
of Quantum Mechanics course). But to get the correct intuition, we can think in 
the following way: in general the solid will have many phonons inside it and each of 
these phonons could sit in any one of the single-phonon states that we counted above. 
Suppose that there are three phonons sitting in the same state. We say that this state 
is occupied three times. In this language, each of the states above can be occupied 
an arbitrary number of times, restricted only by the energy available. If you want to 
describe the total state of the system, you need to say how many phonons are in the 



g{uj)duj 



2n 2 c 2 s 




-70- 



first state, and how many are in the second state and so on. The number of one-phonon 
states is then playing the role of the number of degrees of freedom: it is the number of 
things you can excite. 

The net result of this argument is to equate 

Ar Vu 3 fQn 2 N\ 1/3 

We see that up is related to the atomic spacing (V/N) 1 ^ 3 as we anticipated above, 
but now we have the coefficient too. However, in some sense, the argument of Debye is 
"answer analysis" . It was constructed to ensure that we have the right high-temperature 
behaviour that we will see below. 

From the maximum frequency u D we can construct an associated energy scale, hup, 
and temperature scale, 

rp hujf) 

J-D = -j 

k B 

This is known as the Debye temperature. It provides a way of characterising all solids: 
it is the temperature at which the highest frequency phonon starts to become excited. 
To ranges from around 100K for very soft materials such as lead through to 2000fT 
for hard materials such as diamond. Most materials have Debye temperatures around 
room temperature (±100K or so). 

Heat Capacity of Solids 

All that remains is to put the pieces together. Like photons, the number of phonons is 
not conserved. The partition function for phonons of a fixed frequency, u>, is the same 
as for photons (3.7), 

Summing over all frequencies, the partition function is then 



logZ phonon = / dcu g(uj) log 
Jo 



where the partition function Z u for a single phonon of frequency u is the same as that 
of a photon (3.7). The total energy in sound waves is therefore 

hwaiuj) 3Vh r D , lu 3 
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Figure 17: Experimental data for heat capacities. The solid line is the Debye prediction. 
(Source: D. Schroeder An Introduction to Thermal Physics) 

We again rescale the integration variable to x = /3hw so the upper limit of the 
integral becomes x^ = Tjj/T. Then we have 

QT/ fT D /T 3 

E= n Z ^ M 4 / dx 
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The integral is a function of Tp/T. It has no analytic expression. However, we can 
look in the two extremes. Firstly, for T <C we can replace the upper limit of the 
integral by infinity. We're then left with the same definite integral that we appeared 
for photons, I = f°° dx x 3 / (e x — 1) = 7r 4 /15. In this low-temperature regime, the heat 
capacity is proportional to T 3 , 

c -f^w r3 {t<<Td) (3 ' 13) 

It is often expressed in terms of the Debye temperature Tp, so it reads 

c '="*'Hr©' (3 ' 14) 

In contrast, at temperatures T ^> we only integrate over small values of x, allowing 
us to Taylor expand the integrand, 

rT D /T 3 f T D /T 1 / Tn \3 
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This ensures that the energy grows linearly with T and the heat capacity tends towards 
a constant value 

Vk 4 T 3 

Cv = 7r^TJ = mk B (T»T D ) (3.15) 

This high-temperature behaviour has been known experimentally since the early 1800 's. 
It is called the Dulong-Petit law. Debye's argument for the value of ujo was basically 
constructed to reproduce the coefficient 3N in the formula above. This was known 
experimentally, but also from an earlier model of vibrations in a solid due to Einstein. 
(You met the Einstein model in the first problem sheet). Historically, the real suc- 
cess of the Debye model was the correct prediction of the T 3 behaviour of Cy at low 
temperatures. 

In most materials the heat capacity is dominated by the phonon contribution. (In 
metals there is an additional contribution from conduction electrons that we will cal- 
culate in Section 3.6). The heat capacity of three materials is shown in Figure 17, 
together with the predictions from the Debye model. As you can see, it works very 
well! The deviation from theory at high temperatures is due to differences between Cy 
and C p , the heat capacity at constant pressure. 



What's Wrong with the Debye Model? 

As we've seen, the Debye model is remarkably accurate in capturing the heat capacities 
of solids. Nonetheless, it is a caricature of the physics. The most glaring problem is our 
starting point (3.12). The relationship E = hu between energy and frequency is fine; 
the mistake is the relationship between frequency u and wavevector (or momentum) 
k, namely u = kc s . Equations of this type, relating energy and momentum, are called 
dispersion relations. It turns out that that the dispersion relation for phonons is a little 
more complicated. 



It is not hard to compute the dispersion relation for phonons. 
(You will, in fact, do this calculation in the Applications of 
Quantum Mechanics course). For simplicity, we'll work with X.\' 
a one dimensional periodic lattice of N atoms as shown in 
the figure. The equilibrium position of each atom is x\ = la 
and we impose periodic boundary conditions by insisting that 
xm+i = X\. Let ui be the deviation from equilibrium, ui = 
X\ — la. If we approximate the bonds joining the atoms as har- Figure 18: 



monic oscillators, the Hamiltonian governing the vibrations is 

H = ^ + f ~ Ul+l)2 

i i 

where a is a parameter governing the strength of the bonds between atoms. The 
equation of motion is 

a 

Ui = —{2ui - u i+ i - Ui-i) 
m 

This is easily solved by the discrete Fourier transform. We make the ansatz 
Plugging this into the equation of motion gives the dispersion relation 




To compute the partition function correctly in this model, we would have to revisit the 
density of states using the new dispersion relation E(k) = ftw k . The resulting integrals 
are messy. However, at low temperatures only the smallest frequency modes are excited 
and, for small ka, the sin function is approximately linear. This means that we get 
back to the dispersion relation that we used in the Debye model, oo = kc s , with the 
speed of sound given by c s = a^Ja/m. Moreover, at very high temperatures it is simple 
to check that this model gives the Dulong-Petit law as expected. It deviates from the 
Debye model only at intermediate temperatures and, even here, this deviation is mostly 
negligible. 



3.4 The Diatomic Gas Revisited 

With a bit of quantum experience under our belt, we can look again at the diatomic 
gas that we discussed in Section 2.4. Recall that the classical prediction for the heat 
capacity — Cy = \Nk B — only agrees with experiment at very high temperatures. 
Instead, the data suggests that as the temperature is lowered, the vibrational modes 
and the rotational modes become frozen out. But this is exactly the kind of behaviour 
that we expect for a quantum system where there is a minimum energy necessary to 
excite each degree of freedom. Indeed, this "freezing out" of modes saved us from 
ultra-violet catastrophe in the case of blackbody radiation and gave rise to a reduced 
heat capacity at low temperatures for phonons. 
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Let's start with the rotational modes, described by the Hamiltonian (2.18). Treating 
this as a quantum Hamiltonian, it has energy levels 

E = ^j(j + 1) j=0,l,2,... 

The degeneracy of each energy level is 2j + 1. Thus the rotational partition function 
for a single molecule is 



Z Tot = ^(2j + IK 
3=0 



When T 3> h 2 /2IkB, we can approximate the sum by the integral to get 
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which agrees with our result for the classical partition function (2.19). 

In contrast, for T <C h 2 /2IkB all states apart from j = effectively decouple and 
we have simply Z TOt ~ 1. At these temperatures, the rotational modes are frozen at 
temperatures accessible in experiment so only the translational modes contribute to 
the heat capacity. 

This analysis also explains why there is no rotational contribution to the heat capacity 
of a monatomic gas. One could try to argue this away by saying that atoms are point 
particles and so can't rotate. But this simply isn't true. The correct argument is that 
the moment of inertia I of an atom is very small and the rotational modes are frozen. 
Similar remarks apply to rotation about the symmetry axis of a diatomic molecule. 

The vibrational modes are described by the harmonic oscillator. You already com- 
puted the partition function for this on the first examples sheet (and, in fact, implicitly 
in the photon and phonon calculations above). The energies are 

E = hw(n+ i) 

and the partition function is 

vlb ^ ' ^ 1-e-P** 2sinh(/^/2) 

n n mi/ 

At high temperatures f3tvu) 1, we can approximate the partition function as Z vih m 
1/ (3hw which again agrees with the classical result (2.20). At low temperatures (3fko 3> 
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1, the partition function becomes Z^ e -^/2_ This is a contribution from the zero- 
point energy of the harmonic oscillator. It merely gives the expected additive constant 
to the energy per particle, 

Evih = - log Z vib « — 

and doesn't contribute the heat capacity. Once again, we see how quantum effects 
explain the observed behaviour of the heat capacity of the diatomic gas. The end 
result is a graph that looks like that shown in Figure 11. 

3.5 Bosons 

For the final two topics of this section, we will return again to the simple monatomic 
ideal gas. The classical treatment that we described in Section 2.2 has limitations. As 
the temperature decreases, the thermal de Broglie wavelength, 

V mk sT 

gets larger. Eventually it becomes comparable to the inter-particle separation, (V/N) 1 ^. 
At this point, quantum effects become important. If the particles are non-interacting, 
there is really only one important effect that we need to consider: quantum statistics. 

Recall that in quantum mechanics, particles come in two classes: bosons and fermions. 
Which class a given particle falls into is determined by its spin, courtesy of the spin- 
statistics theorem. Integer spin particles are bosons. This means that any wavefunction 
must be symmetric under the exchange of two particles, 

tp(ri,r 2 ) = ^(f 2 ,fi) 

Particles with —integer spin are fermions. They have an anti-symmetrized wavefunc- 
tion, 

ip(n,r 2 ) = -^(f 2 ,fi) 

At low temperatures, the behaviour of bosons and fermions is very different. All familiar 
fundamental particles such as the electron, proton and neutron are fermions. But an 
atom that contains an even number of fermions acts as a boson as long as we do not 
reach energies large enough to dislodge the constituent particles from their bound state. 
Similarly, an atom consisting of an odd number of electrons, protons and neutrons will 
be a fermion. (In fact, the proton and neutron themselves are not fundamental: they 
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are fermions because they contain three constituent quarks, each of which is a fermion. 
If the laws of physics were different so that four quarks formed a bound state rather 
than three, then both the proton and neutron would be bosons and, as we will see in 
the next two sections, nuclei would not exist!). 

We will begin by describing the properties of bosons and then turn to a discussion 
of fermions to Section 3.6. 

3.5.1 Bose-Einstein Distribution 

We'll change notation slightly from earlier sections and label the single particle quantum 
states of the system by |r). (We used \n) previously, but n will be otherwise occupied 
for most of this section). The single particle energies are then E r and we'll assume that 
our particles are non-interacting. In that case, you might think that to specify the state 
of the whole system, you would need to say which state particle 1 is in, and which state 
particle 2 is in, and so on. But this is actually too much information because particle 1 
and particle 2 are indistinguishable. To specify the state of the whole system, we don't 
need to attach labels to each particle. Instead, it will suffice to say how many particles 
are in state 1 and how many particles are in state 2 and so on. 

We'll denote the number of particles in state |r) as n r . If we choose to work in the 
canonical ensemble, we must compute the partition function, 



where the sum is over all possible ways of partitioning N particles into sets {n r } subject 



makes the sums tricky. This means that the canonical ensemble is rather awkward 
when discussing indistinguishable particles. It turns out to be much easier to work in 
the grand canonical ensemble where we introduce a chemical potential fi and allow the 
total number of particles N to fluctuate. 

Life is simplest it we think of each state |r) in turn. In the grand canonical ensemble, 
a given state can be populated by an arbitrary number of particles. The grand partition 
function for this state is 




{n r } 



to the constraint that J2 r n r = N. Unfortunately, the need to impose this constraint 




1 




Notice that we've implicitly assumed that the sum above converges, which is true only 
if (E r — jj) > 0. But this should be true for all states E r . We will set the ground state 



to have energy E = 0, so the grand partition function for a Bose gas only makes sense 
if 



/i<0 (3.16) 

Now we use the fact that all the occupation of one state is independent of any other. 
The full grand partition function is then 



e -P(E r -p) 



From this we can compute the average number of particles, 
N = \h ^ = ^ eP&-* - 1 " 

r r 

Here (n r ) denotes the average number of particles in the state |r), 

<*> = ^ _ l 

This is the Bose-Einstein distribution. In what follows we will always be interesting 
in the thermodynamic limit where fluctuations around the average are negligible. For 
this reason, we will allow ourselves to be a little sloppy in notation and we write the 
average number of particles in |r) as n r instead of (n r ). 

Notice that we've seen expressions of the form (3.17) already in the calculations for 
photons and phonons — see, for example, equation (3.7). This isn't coincidence: both 
photons and phonons are bosons and the origin of this term in both calculations is the 
same: it arises because we sum over the number of particles in a given state rather 
than summing over the states for a single particle. As we mentioned in the Section 
on blackbody radiation, it is not really correct to think of photons or phonons in the 
grand canonical ensemble because their particle number is not conserved. Nonetheless, 
the equations are formally equivalent if one just sets /i — 0. 

In what follows, it will save ink if we introduce the fugacity 

z = (3.18) 

Since fx < 0, we have < z < 1. 
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Ideal Bose Gas 

Let's look again at a gas of non-relativistic particles, now through the eyes of quantum 
mechanics. The energy is 

2m 

As explained in Section 3.1, we can replace the sum over discrete momenta with an 
integral over energies as long as we correctly account for the density of states. This 
was computed in (3.2) and we reproduce the result below: 

(3 - 19) 

From this, together with the Bose-Einstein distribution, we can easily compute the 
total number of particles in the gas, 



TV 



There is an obvious, but important, point to be made about this equation. If we can do 
the integral (and we will shortly) we will have an expression for the number of particles 
in terms of the chemical potential and temperature: N = N(fi,T). That means that if 
we keep the chemical potential fixed and vary the temperature, then N will change. But 
in most experimental situations, N is fixed and we're working in the grand canonical 
ensemble only because it is mathematically simpler. But nothing is for free. The price 
we pay is that we will have to invert equation (3.20) to get it in the form \x = fJ,(N, T). 
Then when we change T, keeping N fixed, /i will change too. We have already seen an 
example of this in the ideal gas where the chemical potential is given by (2.14). 

The average energy of the Bose gas is, 

And, finally, we can compute the pressure. In the grand canonical ensemble this is, 

P V = ^\ogZ = -^J dEg(E)log(l-ze-P E ) 

We can manipulate this last expression using an integration by parts. Because g(E) ~ 
E 1 / 2 , this becomes 



P V 
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This is implicitly the equation of state. But we still have a bit of work to do. Equation 
(3.21) gives us the energy as a function of \x and T. And, by the time we have inverted 
(3.20), it gives us /i as a function of N and T. Substituting both of these into (3.22) 
will give the equation of state. We just need to do the integrals. . . 



3.5.2 A High Temperature Quantum Gas is (Almost) Classical 

Unfortunately, the integrals (3.20) and (3.21) look pretty fierce. Shortly, we will start 
to understand some of their properties, but first we look at a particularly simple limit. 
We will expand the integrals (3.20), (3.21) and (3.22) in the limit z = < 1. We'll 
figure out the meaning of this expansion shortly (although there's a clue in the title of 
this section if you're impatient). Let's look at the particle density (3.20), 

N 1 (2m\ 3/2 f°° , E 1 / 2 





(2m 
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4vr 2 1 





V 4tt 2 V h 2 J J z~ l e^ E - 1 
, 3/2 r lrn ze-^E 1 / 2 





)3/2 «oo 
J dx ^e- x (l + ze- x + ...) 

where we made the simple substitution x = f3E. The integrals are all of the Gaussian 
type and can be easily evaluated by making one further substitution x = u 2 . The final 
answer can be conveniently expressed in terms of the thermal wavelength A, 

Now we've got the answer, we can ask what we've done! What kind of expansion is z 
1? From the above expression, we see that the expansion is consistent only if X 3 N/V <sC 
1, which means that the thermal wavelength is much less than the interparticle spacing. 
But this is true at high temperatures: the expansion 2 < 1 is a high temperature 
expansion. 



At first glance, it is surprising that z = e^ 4 1 corresponds to high temperatures. 
When T — > 00, (3 — > so naively it looks as if z — > 1. But this is too naive. If we keep 
the particle number N fixed in (3.20) then \x must vary as we change the temperature, 
and it turns out that fi — > — 00 faster than (5 — > 0. To see this, notice that to leading 
order we need z/X 3 to be constant, soz~ T~ 3//2 . 
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High Temperature Equation of State of a Bose Gas 

We now wish to compute the equation of state. We know from (3.22) that pV = ^E. 
So we need to compute E using the same z expansion as above. From (3.21), the 
energy density is 

E 1 (2m\ 3/2 f°° , p £ 3/2 

3Z (3.24) 



2A 3 /3 V 4^ 

The next part's a little fiddly. We want to eliminate z from the expression above in 
favour of N/V. To do this we invert (3.23), remembering that we're working in the 
limit z < 1 and X 3 N/V < 1. This gives 

X 3 N ( 1 X 3 N 

1 = —— + ... 



V V 2^ V 
which we then substitute into (3.24) to get 

^ 3N f 1 X 3 N \ ( 1 X 3 N 

and finally we substitute this into pV = |E to get the equation of state of an ideal 
Bose gas at high temperatures, 

which reproduces the classical ideal gas, together with a term that we can identify 
as the second virial coefficient in the expansion (2.21). However, this correction to 
the pressure hasn't arisen from any interactions among the atoms: it is solely due to 
quantum statistics. We see that the effect of bosonic statistics on the high temperature 
gas is to reduce the pressure. 

3.5.3 Bose-Einstein Condensation 

We now turn to the more interesting question: what happens to a Bose gas at low 
temperatures? Recall that for convergence of the grand partition function we require 
/i < so that z = G (0,1). Since the high temperature limit is z — > 0, we'll 
anticipate that quantum effects and low temperatures come about as z — > 1. 
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Recall that the number density is given by (3.20) which we write as 
N 1 (2mk B T\* 12 f°° , x 1 ' 2 1 



V 4tt 2 



( 2mk B T \^ [°° x^ 2 1 



where the function (73/2(2) is one of a class of functions which appear in the story of 
Bose gases, 



x n ~ 1 



9n(z) = r—— / dx (3.27) 

Y[n) J z 1 e x - 1 

These functions are also known as polylogarithms and sometimes denoted as Li n (z) = 
g n (z). The function g 3 / 2 is relevant for the particle number and the function g 5 / 2 
appears in the calculation of the energy in (3.21). For reference, the gamma function 
has value T(3/2) = y/n/2. (Do not confuse these functions g n (z) with the density of 
states g(E). They are not related; they just share a similar name). 

In Section 3.5.2 we looked at (3.26) in the T — > 00 limit. There we saw that A — > 
but the function #3/2(2) — > in just the right way to compensate and keep N/V fixed. 
What now happens in the other limit, T — > and A — > 00. One might harbour the hope 
that #3/2(2) — > 00 and again everything balances nicely with N/V remaining constant. 
We will now show that this can't happen. 

A couple of manipulations reveals that the integral (3.27) can be expressed in terms 
of a sum, 
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But the integral that appears in the last line above is nothing but the definition of the 
gamma function T(n). This means that we can write 



x z m 



9»W = T,—n ( 3 - 28 ) 



m 
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We see that g n (z) is a monotonically increasing function of z. Moreover, at z — 1, it is 
equal to the Riemann zeta function 

9n{±) = C(n) 

For our particular problem it will be useful to know that C(3/2) ~ 2.612. 

Let's now return to our story. As we decrease T in (3.26), keeping N/V fixed, z and 
hence g 3 / 2 (z) must both increase. But z can't take values greater than 1. When we 
reach z — 1, let's denote the temperature as T = T c . We can determine T c by setting 
z — 1 in equation (3.26), 

T ^(^)(l"Y 3 (3 , 9) 



\k B m J V.C(3/2) V 

What happens if we try to decrease the temperature below T c ? Taken at face value, 
equation (3.26) tells us that the number of particles should decrease. But that's a 
stupid conclusion! It says that we can make particles disappear by making them cold. 
Where did they go? We may be working in the grand canonical ensemble, but in 
the thermodynamic limit we expect AiV ~ l/y/N which is too small to explain our 
missing particles. Where did we misplace them? It looks as if we made a mistake in 
the calculation. 

In fact, we did make a mistake in the calculation. It happened right at the beginning. 
Recall that back in Section 3.1 we replaced the sum over states with an integral over 
energies, 



Because of the weight V E in the integrand, the ground state with E — doesn't 
contribute to the integral. We can use the Bose-Einstein distribution (3.17) to compute 
the number of states that we expect to be missing because they sit in this E = ground 

state, 

n = -z^—r (3-30) 
z 1 — 1 

For most values of z G (0, 1), there are just a handful of particles sitting in this lowest 
state and it doesn't matter if we miss them. But as z gets very very close to 1 (meaning 
z ~ 1 — 1/ N) then we get a macroscopic number of particles occupying the ground state. 
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It is a simple matter to redo our calculation taking into account the particles in the 
ground state. Equation (3.26) is replaced by 

Now there's no problem keeping N fixed as we take z close to 1 because the additional 
term diverges. This means that if we have finite N, then as we decrease T we can never 
get to z — 1. Instead, z must level out around z « 1 — 1/N as T — > 0. 

For T < T c , the number of particles sitting in the ground state is 0(M). Some simple 
algebra allows us to determine that fraction of particles in the ground state is 

S-X-JLfP/^l-g)'" (3,1) 

At temperatures T < T c , a macroscopic number of atoms discard their individual 
identities and merge into a communist collective, a single quantum state so large that 
it can be seen with the naked eye. This is known as the Bose-Einstein condensate. It 
provides an exquisitely precise playground in which many quantum phenomena can be 
tested. 

Bose-Einstein Condensates in the Lab 

Bose-Einstein condensates (often shortened to BECs) of 
weakly interacting atoms were finally created in 1995, 
some 70 years after they were first predicted. These first 
BECs were formed of Rubidium, Sodium or Lithium 
and contained between N ~ 10 4 — > 10 7 atoms. The Figure 19: UFO=BEC 
transition temperatures needed to create these conden- 
sates are extraordinarily small, around T c ~ 10~ 7 K. Figure 20 shows the iconic colour 
enhanced plots that reveal the existence of the condensate. To create these plots, the 
atoms are stored in a magnetic trap which is then turned off. A picture is taken of the 
atom cloud a short time t later when the atoms have travelled a distance Hkt/m. The 
grey UFO-like smudges above are the original pictures. From the spread of atoms, the 
momentum distribution of the cloud is inferred and this is what is shown in Figure 20. 
The peak that appears in the last two plots reveals that a substantial number of atoms 
were indeed sitting in the momentum ground state. (This is not a k = state because 
of the finite trap and the Heisenberg uncertainty relation). The initial discoverers of 
BECs, Eric Cornell and Carl Wieman from Boulder and Wolfgang Ketterle from MIT, 
were awarded the 2001 Nobel prize in physics. 
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Figure 20: The velocity distribution of Rubidium atoms, taken from the Ketterle lab at 
MIT. The left-hand picture shows T > T c , just before the condensate forms. The middle and 
right-hand pictures both reveal the existence of the condensate. 



Low Temperature Equation of State of a Bose Gas 

The pressure of the ideal Bose gas was computed in (3.22). We can express this in 
terms of our new favourite functions (3.27) as 

2E k B T . . . 
P= 3^7= -^-95/2{z) (3.32) 

Formally there is also a contribution from the ground state, but it is log(l — z)/V 
which is a factor of N smaller than the term above and can be safely ignored. At low 
temperatures, T < T c , we have z ~ 1 and 

P = ^C(5/2) 

So at low temperatures, the equation of state of the ideal Bose gas is very different from 
the classical, high temperature, behaviour. The pressure scales as p ~ T 5 / 2 (recall that 
there is a factor of T 3 / 2 lurking in the A) . More surprisingly, the pressure is independent 
of the density of particles N/V. 

3.5.4 Heat Capacity: Our First Look at a Phase Transition 

Let's try to understand in more detail what happens as we pass through the critical 
temperature T c . We will focus on how the heat capacity behaves on either side of the 
critical temperature. 
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We've already seen in (3.32) that we can express the energy in terms of the function 

#5/2(2), 

E 3k B T 



V 2 A 3 

so the heat capacity becomes 



■#3/2(2) 



Cy _ 1 dE _ 15k B / \ 3 fc B T d# 5/2 rfz 
V~ " V dT ~ -&F9W Z > + 2— — dT (3 - 33) 
The first term gives a contribution both for T < T c and for T > T c . However, the second 
term includes a factor of dz/dT and z is a very peculiar function of temperature: for 
T > T c , it is fairly smooth, dropping off at T = T c . However, as T — > T c , the fugacity 
rapidly levels off to the value z ~ 1 — 1/N. For T < T c , z doesn't change very much at 
all. The net result of this is that the second term only contributes when T > T c . Our 
goal here is to understand how this contribution behaves as we approach the critical 
temperature. 

Let's begin with the easy bit. Below the critical temperature, T < T c , only the first 
term in (3.33) contributes and we may happily set z — 1. This gives the heat capacity, 

1 

Cv = ^yC(5/2) ~ T 3 / 2 (3.34) 

Now we turn to T > T c . Here we have z < 1, so ^5/2(2) < #5/2(1). We also have 
dz/dT < 0. This means that the heat capacity decreases for T > T c . But we know 
that for T < T c , Cy ~ T 3 / 2 so the heat capacity must have a maximum at T = T c . 
Our goal in this section is to understand a little better what the function Cy looks like 
in this region. 

To compute the second term in (3.33) we need to understand both g^ 2 and how z 
changes with T as we approach T c from above. The first calculation is easy is we use 
our expression (3.28), 

9n{z) = dz~ 9n ^ = ~9n-i{z) (3-35) 



m=l 



As T — > T c from above, dg^/2/dT — > ^(3/2), a constant. All the subtleties lie in the 
remaining term, dz/dT. After all, this is the quantity which is effectively vanishing for 
T < T c . What's it doing at T > T c ? To figure this out is a little more involved. We 
start with our expression (3.26), 

7VA 3 

93/2(z) = — T>T C (3.36) 
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and we'll ask what happens to the function #3/2(2) as z — > 1, keeping N fixed. We 
know that exactly at z — 1, #3/2(1) = £(3/2). But how does it approach this value? To 
answer this, it is actually simplest to look at the derivative dg 3 / 2 /dz = gi/2/z, where 



9l/2{z) = Wml dx 



x- 1 ' 2 

z -l e x _ I 



The reason for doing this is that g±/2 diverges as z — > 1 and it is generally easier to 
isolate divergent parts of a function than some finite piece. Indeed, we can do this 
straightforwardly for g±/ 2 by looking at the integral very close to x — 0, we can write 

1 r x- 1 / 2 

9i/2iz) = rm L dx z -\i +X )-i + finite 

dX r h • • • 



r(i/2)7 " (i-z) + x 
2z 1 r , 1 

du + . . . 



^r~z r(i/2) J i + u 



where, in the last line, we made the substution u = \/x/(l — z). So we learn that 
as z — > 1, gi/2(z) — > z{\ — z)~ x l 2 . But this is enough information to tell us how #3/2 
approaches its value at z — 1: it must be 

03/2OO*C(3/2)+A(l-*) 1/2 + ... 

for some constant A. Inserting this into our equation (3.36) and rearranging, we find 
that as T — > T c from above, 

T-T c 



\ Vrj 

1 - 5 ' ' 



T, 



c 



where, in the second line, we used the expression of the critical temperature (3.29). 
B is some constant that we could figure out with a little more effort, but it won't be 
important for our story. From the expression above, we can now determine dz/dT as 
T — > T c . We see that it vanishes linearly at T = T c . 
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Putting all this together, we can determine the expression for the heat capacity 
(3.33) when T > T c . We're not interested in the coefficients, so we'll package a bunch 
of numbers of order 1 into a constant b and the end result is 



C 



v 



15Vk B (T 
-J^9,/2{z) - b 



T, 



The first term above goes smoothly over to the 
expression (3.34) for Cy when T < T c . But 
the second term is only present for T > T c . 
Notice that it goes to zero as T — > T c , which 
ensures that the heat capacity is continuous at 
this point. But the derivative is not continu- 
ous. A sketch of the heat capacity is shown in 
the figure. 




3Nk B /2 



Figure 21: Heat Capacity for a BEC 



Functions in physics are usually nice and 
smooth. How did we end up with a discontinuity in the derivative? In fact, if we 
work at finite N, strictly speaking everything is nice and smooth. There is a similar 
contribution to dz/dT even at T < T c . We can see that by looking again at the 
expressions (3.30) and (3.31), which tell us 



3=1 + 



= 1 + 



1 



(T < T c ) 



n / \ Nl-{T/T c fl\ 

The difference is that while dz/dT is of order one above T c , it is of order 1/N below 
T c . In the thermodynamic limit, iV — > oo, this results in the discontinuity that we saw 
above. This is a general lesson: phase transitions with their associated discontinuities 
can only arise in strictly infinite systems. There are no phase transitions in finite 
systems. 

Superfluid Helium-4 

A similar, but much more pronounced, discontinuity is seen in Helium-4 as it 
becomes a superfluid, a transition which occurs at 2.17 K. The atom contains two 
protons, two neutrons and two electrons and is therefore a boson. (In contrast, Helium- 
3 contains just a single neutron and is a fermion). The experimental data for the heat 
capacity of Helium-4 is shown on the right. The successive graphs are zooming in 
on the phase transistion: the scales are (from left to right) Kelvin, milliKelvin and 
microKelvin. The discontinuity is often called the lambda transition on account of the 
shape of this graph. 
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There is a close connection between Bose- 
Einstein condensation described above and su- 
perfluids: strictly speaking a non-interacting Bose- 
Einstein condensate is not a superfluid but su- 
perfluidity is a consequence of arbitrarily weak 
repulsive interactions between the atoms. How- 
ever, in He-4, the interactions between atoms 




are strong and the system cannot be described 
using the simple techniques developed above. 



Figure 22: 4 He. 



Something very similar to Bose condensation 
also occurs in superconductivity and superfluidity of Helium-3. Now the primary char- 
acters are fermions rather than bosons (electrons in the case of superconductivity). As 
we will see in the next section, fermions cannot condense. But they may form bound 
states due to interactions and these effective bosons can then undergo condensation. 

3.6 Fermions 

For our final topic, we will discuss the fermion gases. Our analysis will focus solely 
on non-interacting fermions. Yet this simple model provides a (surprisingly) good 
first approximation to a wide range of systems, including electrons in metals at low 
temperatures, liquid Helium-3 and white dwarfs and neutron stars. 

Fermions are particles with —integer spin. By the spin-statistics theorem, the wave- 
function of the system is required to pick up a minus sign under exchange of any 
particle, 



As a corollory, the wavefunction vanishes if you attempt to put two identical fermions 
in the same place. This is a reflection of the Pauli exclusion principle which states that 
fermions cannot sit in the same state. We will see that the low-energy physics of a gas 
of fermions is entirely dominated by the exclusion principle. 

We work again in the grand canonical ensemble. The grand partition function for a 
single state |r) is very easy: the state is either occupied or it is not. There is no other 
option. 



V>(fi,f 2 ) 
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So, the grand partition function for all states is Z = f| r Z r , from which we can compute 
the average number of particles in the system 

N = S e /3(E r -^) + I =^2 n r 



r 



where the average number of particles in the state |r) is 

This is the Fermi-Dirac distribution. It differs from the Bose-Einstein distribution only 
by the sign in the denominator. Note however that we had no convergence issues in 
defining the partition function. Correspondingly, the chemical potential /x can be either 
positive or negative for fermions. 

3.6.1 Ideal Fermi Gas 

We'll look again at non-interacting, non-relativistic particles with E = h 2 k 2 /2m . Since 
fermions necessarily have —integer spin, s, there is always a degeneracy factor when 
counting the number of states given by 

9s = 2s + 1 

For example, electrons have spin \ and, correspondingly have a degeneracy of g s = 2 
which just accounts for "spin up" and "spin down" states. We saw similar degeneracy 
factors when computing the density of states for photons (which have two polarizations) 
and phonons (which had three). For non-relativistic fermions, the density of states is 

We'll again use the notation of fugacity, z = e 13 ^. The particle number is 

N = J dE "F^TI (3 - 38) 



The average energy is 



z-^ePz + 1 
And the pressure is 

P V = ^JdE g{E) log (1 + ze-* E ) = ^E (3.39) 
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At high temperatures, it is simple to repeat the steps of Section 3.5.2. (This is one of 
the questions on the problem sheet). Ony a few minus signs differ along the way and 
one again finds that for z <C 1, the equation of state reduces to that of a classical gas, 



pV = Nk B T 1 + 



A 3 iV 



(3.40) 



4y/2g a V 



+ 



Notice that the minus signs filter down to the final answer: the first quantum correction 
to a Fermi gas increases the pressure. 

3.6.2 Degenerate Fermi Gas and the Fermi Surface 

In the extreme limit T — > 0, the Fermi-Dirac distribution becomes very simple: a state 
is either filled or empty. 



It's simple to see what's going on here. Each fermion that we throw into the system 
settles into the lowest available energy state. These are successively filled until we run 
out of particles. The energy of the last filled state is called the Fermi energy and is 
denoted as Ep. Mathematically, it is the value of the chemical potential at T = 0, 



Filling up energy states with fermions is just like throwing balls into a box. With one 
exception: the energy states of free particles are not localised in position space; they 
are localised in momentum space. This means that successive fermions sit in states 
with ever-increasing momentum. In this way, the fermions fill out a ball in momentum 
space. The momentum of the final fermion is called the Fermi momentum and is related 
to the Fermi energy in the usual way: fikp = (2mEp) l l 2 . All states with wavevector 
\k\ < kp are filled and are said to form the Fermi sea or Fermi sphere. Those states 
with \k\ = kp lie on the edge of the Fermi sea. They are said to form the Fermi 
surface. The concept of a Fermi surface is extremely important in later applications to 
condensed matter physics. 

We can derive an expression for the Fermi energy in terms of the number of particles 
N in the system. To do this, we should appreciate that we've actually indulged in a 
slight abuse of notation when writing (3.41). In the grand canonical ensemble, T and 
fi are independent variables: they're not functions of each other! What this equation 
really means is that if we want to keep the average particle number N in the system 




for E < /j, 
for E > fi 



»(T = 0) = Ep 



(3.41) 
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fixed (which we do) then as we vary T we will have to vary fi to compensate. So a 
slightly clearer way of defining the Fermi energy is to write it directly in terms of the 
particle number 



"=r<W = £f¥fV (3.42) 



Or, inverting, 



EF *(<?l»f (3.43) 
2m V 9s V) y ! 



The Fermi energy sets the energy scale for the system. There is an equivalent tem- 
perature scale, Tp = Ep/ks- The high temperature expansion that resulted in the 
equation of state (3.40) is valid at temperatures T > T F . In contrast, temperatures 
T < T F are considered "low" temperatures for systems of fermions. Typically, these 
low temperatures do not have to be too low: for electrons in a metal, T F ~ 10 4 fT; for 
electrons in a white dwarf, T F > 10 7 K. 

While E F is the energy of the last occupied state, the average energy of the system 
can be easily calculated. It is 

i-E F o 

E = / dE Eg(E) = -NE F (3.44) 
Jo 5 

Similarly, the pressure of the degenerate Fermi gas can be computed using (3.39), 

pV = -NE F (3.45) 
5 

Even at zero temperature, the gas has non-zero pressure, known as degeneracy pressure. 
It is a consequence of the Pauli exclusion principle and is important in the astrophysics 
of white dwarf stars and neutron stars. (We will describe this application in Section 
3.6.5). The existence of this residual pressure at T = is in stark contrast to both the 
classical ideal gas (which, admittedly, isn't valid at zero temperature) and the bosonic 
quantum gas. 

3.6.3 The Fermi Gas at Low Temperature 

We now turn to the low-temperature behaviour of the Fermi gas. As mentioned above, 
"low" means T <ti Tp, which needn't be particularly low in everyday terms. The 
number of particles N and the average energy E are given by, 

9(E) 



DO 



N = I dE ^ (3.46) 



o 
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Figure 23: The Fermi-Dirac Distribution function at T = and very small T. The distri- 
bution differs from the T = ground state only for a range of energies ksT around Ep. 



Our goal is to firstly understand how the chemical potential //, or equivalently the 
fugacity z = e 13 ^, changes with temperature when N is held fixed. From this we can 
determine how E changes with temperature when N is held fixed. 

There are two ways to proceed. The first is a direct approach on the problem by 
Taylor expanding (3.46) and (3.47) for small T. But it turns out that this is a little 
delicate because the starting point at T = involves the singular distribution shown 
on the left-hand side of Figure 23 and it's easy for the physics to get lost in a morass of 
integrals. For this reason, we start by developing a heuristic understanding of how the 
integrals (3.46) and (3.47) behave at low temperatures which will be enough to derive 
the required results. We will then give the more rigorous expansion - sometimes called 
the Sommerfeld expansion - and confirm that our simpler derivation is indeed correct. 

The Fermi-Dirac distribution (3.37) at small temperatures is sketched on the right- 
hand side of Figure 23. The key point is that only states with energy within ksT of 
the Fermi surface are affected by the temperature. We want to insist that as we vary 
the temperature, the number of particles stays fixed which means that dN/dT = 0. 
I claim that this holds if, to leading order, the chemical potential is independent of 
temperature, so 



and 




(3.47) 



where, for non-relativistic fermions, the density of states is 




d/j, 







(3.48) 



dT 



T=0 
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Let's see why this is the case. The change in particle number can be written as 

dN r dE 9(E) 



dT dT J e^ E ~ri + 1 
, " 00 d ( 1 

,„ dE ' AE) w {^wn 



/ 

Jo 

g(E F ) [ 

Jo 



dT \ e (KE-E F ) + l 



There are two things going on in the step from the second line to the third. Firstly, 
we are making use of the fact that, for ksT <C Ep, the Fermi-Dirac distribution 
only changes significantly in the vicinity of Ep as shown in the right-hand side of 
Figure 23. This means that the integral in the middle equation above is only receiving 
contributions in the vicinity of Ep and we have used this fact to approximate the density 
of states g(E) with its value at g(Ep). Secondly, we have used our claimed result (3.48) 
to replace the total derivative d/dT (which acts on the chemical potential) with the 
partial derivative d/dT (which doesn't) and \i is replaced with its zero temperature 
value E F . 

Explicitly differentiating the Fermi-Dirac distribution in the final line, we have 

This integral vanishes because (E — E F ) is odd around E F while the cosh function 
is even. (And, as already mentioned, the integral only receives contributions in the 
vicinity of Ep). 

Let's now move on to compute the change in energy with temperature. This, of 
course, is the heat capacity. Employing the same approximations that we made above, 
we can write 



BE 

Cv ~df 



N,V 
oo 

dE 



r°° d 



o 







Epg(E F ) + ^g(Ep)(E-E f ) 



dT V e^ E ~ E E) + 1 



However, this time we have not just replaced Eg(E) by Epg(Ep), but instead Taylor 
expanded to include a term linear in (E — Ep). (The factor of 3/2 comes about because 
Eg(E) ~ E z l 2 . The first E F g(E F ) term in the square brackets vanishes by the same 
even/odd argument that we made before. But the (E — E F ) term survives. 
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Writing x — (5{E — Ep), the integral becomes 

where we've extended the range of the integral from — oo to +00, safe in the knowledge 
that only the region near Ep (or x — 0) contributes anyway. More importantly, how- 
ever, this integral only gives an overall coefficient which we won't keep track of. The 
final result for the heat capacity is 

C v ~Tg(E F ) 

There is a simple way to intuitively understand this linear behaviour. At low temper- 
atures, only fermions within k B T of Ep are participating in the physics. The number 
of these particles is ~ g(E F )k B T. If each picks up energy ~ k B T then the total energy 
of the system scales as E ~ g(E F )(k B T) 2 , resulting in the linear heat capacity found 
above. 

Finally, the heat capacity is often re-expressed in a slightly different form. Using 
(3.42) we learn that iV ~ E 3 J 2 which allows us to write, 

C v ~ Nk B (Jr^J (3-49) 

Heat Capacity of Metals 

One very important place that we can apply the theory above it to metals. We can try 
to view the conduction electrons — those which are free to move through the lattice 
- as an ideal gas. At first glance, this seems very unlikely to work. This is because we 
are neglecting the Coulomb interaction between electrons. Nonetheless, the ideal gas 
approximation for electrons turns out to work remarkably well. 

From what we've learned in this section, we would expect two contributions to the 
heat capacity of a metal. The vibrations of the lattice gives a phonon contribution 
which goes as T 3 (see equation (3.13)). If the conduction electrons act as an ideal gas, 
they should give a linear term. The low-temperature heat capacity can therefore be 
written as 

C v = 7 T + aT 3 

Experimental data is usually plotted as C v /T vs T 2 since this gives a straight line 
which looks nice. The intercept with the Cv axis tells us the electron contribution. 
The heat capacity for copper is shown in the figure. 
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We can get a rough idea of when the phonon and 
electron contributions to heat capacity are compara- 
ble. Equating (3.14) with (3.56), and writing 247r 2 /5 ~ 




50, we have the the two contributions are equal when 
T 2 ~ Tg/501>. Ballpark figures are T D ~ 10 2 K and 
Tp ~ 10 4 if which tells us that we have to get down 
to temperatures of around IK or so to see the electron 
contribution. 



For many metals, the coefficient 7 of the linear heat 
capacity is fairly close to the ideal gas value (within, 
say, 20% or so). Yet the electrons in a metal are far 



Figure 24: The heat capacity 
of copper (taken from A. Tari 
" The Specific Heat of Matter at 
Low Temperatures") 



from free: the Coulomb energy from their interations 
is at least as important as the kinetic energy. So why 

does the ideal gas approximation work so well? This was explained in the 1950s by 
Landau and the resulting theory — usually called Landau's Fermi liquid theory — is 
the basis for our current understanding of electron systems. 

3.6.4 A More Rigorous Approach: The Sommerfeld Expansion 

The discussion in Section 3.6.3 uses some approximations but highlights the physics 
of the low-temperature fermi gas. For completeness, here we present a more rigorous 
derivation of the low-temperature expansion. 



To start, we introduce some new notation for the particle number and energy inte- 
grals, (3.46) and (3.47). We write, 




7§ h/2{z) 



(3.50) 



and 




where A = \j2irh 2 /mkBT is the familiar thermal wavelength and the functions f n (z) 
are the fermionic analogs of g n defined in (3.27), 




(3.52) 
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where it is useful to recall T(3/2) = and T(5/2) = 3^/7^/4 (which follows, of 

course, from the property T(n + 1) = nT(n)). This function is an example of a poly- 
logarithm and is sometimes written as Li n (—z) = —f n (z). 

Expanding f n (z) 

We will now derive the large z expansion of f n (z), sometimes called the Sommerfeld 
expansion. The derivation is a little long. We begin by splitting the dx integral into 
two parts, 



rPii X n ~ f°° 

r(n)f n (z) = dx — j— — — + / dx 
Jo z l e x + l J, 



to z' 1 ^ + 1 

Pn / i \ roc n-1 

dx x n ~ x 1 + / dx — 

(log*)™ [to x n ~ l [°° x n - 1 



n 



o 1 + ze x J/3n ' z 1 e x + l 



We now make a simple change of variable, r/i = (3\i — x for the first integral and 
r/2 = x — f3 ji for the second, 

(log*)* re» (p^-vir- 1 r, (PV+V2T- 1 

r{n)f n (z) = / di] x — — — + / d7] 2 



n Jo 1 + e" 1 Jo 1 + e m 

So far we have only massaged the integral into a useful form. We now make use of the 
approximation (3fi ^> 1. Firstly, we note that the first integral has a e m suppression 
in the denominator. If we replace the upper limit of the integral by 00, the expression 
changes by a term of order e _/3/i = z^ 1 which is a term we're willing to neglect. With 
this approximation, we can combine the two integrals into one, 

(log z) n [°° {(3n + r?)"- 1 - ((3/i - riY 1 - 1 



1 + & 



n Jo 

We now Taylor expand the numerator, 

tfv + vT' 1 - tfv - vY 1 - 1 = (/^r 1 [(i + vl^Y' 1 - (i - v/M 71 - 1 } 

= i^Y' 1 

= 2(ra- l)((3fiY' 2 V + ■■■ 

From which we get 

(logz) B , n ,_ lWl v„- 2 r,„ V 



l + (n-l)^ + ...-l + (n-l)^ + 



T(n)f n (z) = + 2(n - l)(log*)- 2 / 

n Jo 



drj 



e" 
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We're left with a fairly simple integral. 



3 n r ye-v 

dr] = / di] 

en + 1 L 1 + e-v 



poo 00 

/ dwJZ 

J o m =i 



~ (_l)m+l 

' m 2 

m=l 



e —mri /_j\m+l 



ofai we" 



The integral in the last line is very simple: f^°du ue u — 1. We have only to do the 
sum. There are a number of ways to manipulate this. We chose to write 

111 / 111 \ / 1 1 1 

1 ~2 2 " + 3 2 "~42'"~ V 1 + 22 + 32 + 42 + "7 _2 l v 2 2 " + 42 + 6 2 " + "' 

1 -f)-( 1 4 + f + ---) <3 - 53) 

This final sum is well known. It is 1/m 2 = ((2) = 7r 2 /6. (The original proof of this, 
by Euler, proceeds by Taylor expanding sinx/x, writing resulting sum as a product of 
roots, and then equating the x 2 term). 

After all this work, we finally have the result we want. The low temperature expan- 
sion of f n (z) is an expansion in l/logz = 1 //?//, 

,„ W (l + ^rf^ + . . .) (3.54) 

Y[n + 1) V 6 (log^) 2 J 

Back to Physics: Heat Capacity of a Fermi Gas Again 

The expansion (3.54) is all we need to determine the leading order effects of temperature 
on the Fermi gas. The number of particles (3.50) can be expressed in terms of the 
chemical potential as 

N 9 s /0 ^3/2 



V 6ir 2 h 3 




(2m/i)^ N + _ +... (3.55) 



But the number of particles determines the Fermi energy through (3.43). For fixed 
particle number N, we can therefore express the chemical potential at finite temperature 
in terms of the Fermi energy. After a small amount of algebra, we find 
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We see that the chemical potential is a maximum at T = and decreases at higher 
T. This is to be expected: recall that by the time we get to the classical gas, the 
chemical potential must be negative. Moreover, note that the leading order correction 
is quadratic in temperature rather than linear, in agreement with our earlier claim 
(3.48). 



The energy density can be computed from (3.51) and (3.54). It is 



Our real interest is in the heat capacity. However, with fixed particle number N, the 
chemical potential also varies with temperature at quadratic order. We can solve this 
problem by dividing by (3.55) to get 

E IE, f 1 | 5tt 2 fk B T\ 2 



N 5 1 12 V E 

Now the only thing that depends on temperature on the right-hand side is T itself. 
From this we learn that the heat capacity of a low temperature Fermi gas is linear in 
T, 

dE 



Cv ~dT 



7T 2 T 

= Nk B -— (3.56) 



But we now have the correct coefficient that completes our earlier result (3.49). 
3.6.5 White Dwarfs and the Chandrasekhar limit 

When stars exhaust their fuel, the temperature T — > and they have to rely on the Pauli 
exclusion principle to support themselves through the degeneracy pressure (3.45). Such 
stars supported by electron degeneracy pressure are called white dwarfs. In addition 
to the kinetic energy -Ekinctk of fermions given in (3.43) the system has gravitational 
energy. If we make the further assumption that the density in a star of radius R is 
uniform, then 

3 G N M 2 
grav ~ 5 R 

where Gjy is Newton's constant. In the problem set, you will minimise -E^v + -^kinetic 
to find the relationship between the radius and mass of the star, 

R ~ M- 1/3 

This is unusual. Normally if you throw some stuff on a pile, the pile gets bigger. Not 
so for stars supported by degeneracy pressure. The increased gravitational attraction 
wins and causes them to shrink. 
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As the star shrinks, the Fermi energy (3.43) grows. Eventually it becomes comparable 
to the electron mass m e and our non-relativistic approximation breaks down. We can re- 
do our calculations using the relativistic formula for the density of states (3.4) (which 
needs to be multiplied by a factor of 2 to account for the spin degeneracy). In the 
opposite regime of ultra- relativistic electrons, where E ^> m e , then we can expand the 
density of states as 



9(E) 



V 



1T 2 h 3 C 3 



E — 



1 4 

m c 



+ 



from which we can determine the kinetic energy due to fermions, replacing the non- 
relativistic result (3.44) by 



E]<\ 



V 



kinetic 



TT 2 h 3 C 3 



1 4 m 2 c 4 2 
-E F - —E F + ... 



The Fermi energy can be expressed in terms of the particle number by 

V i I 



N 



n 2 h 3 c 3 V 3 



-Ei 



1 4 

m c 



-E F + ... 



The total energy of the system is then given by, 



E. 



grav 



3hc /9ttM 4 ^ 1/3 
~T V 4m 4 



- ^G N M 2 



R 



+ 0(R) 



where m p is the proton mass (and we're taking M = Nm p as a good approximation 
to the full mass of the star). If the term above is positive then the star once again 
settles into a minimum energy state, where the 1/R term balances the term that grows 
linearly in R. But if the 1/R term is negative, we have a problem: the star is unstable 
to gravitational collapse and will shrink to smaller and smaller values of R. This occurs 
when the mass exceeds the Chandrasekhar limit, M > Mq. Neglecting factors of 2 and 
7r, this limit is 



M, 



c 



he \ 



3/2 



Gn J 



mi 



(3.57) 



This is roughly 1.5 times the mass of the Sun. Stars that exceed this bound will not 
end their life as white dwarfs. Their collapse may ultimately be stopped by degeneracy 
pressure of neutrons so that they form a neutron star. Or it may not be stopped at all 
in which case they form a black hole. 
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The factor in brackets in (3.57) is an interesting mix of fundamental constants as- 
sociated to quantum theory, relativity and gravity. In this combination, they define a 
mass scale called the Planck mass, 



When energy densities reach the Planck scale, quantum gravity effects are important 
meaning that we have to think of quantising spacetime itself. We haven't quite reached 
this limit in white dwarf stars. (We're shy by about 25 orders of magnitude or so!). 
Nonetheless the presence of the Planck scale in the Chandra limit is telling us that 
both quantum effects and gravitational effects are needed to understand white dwarf 
stars: the matter is quantum; the gravity is classical. 

3.6.6 Pauli Paramagnetism 

"With heavy heart, I have decided that Fermi-Dirac, not Einstein, is the 
correct statistics, and I have decided to write a short note on paramag- 
netism" 

Wolfgang Pauli, in a letter to Schrodinger, reluctantly admitting that elec- 
trons are fermions 

We now describe a gas of electrons in a background magnetic field B. There are 
two effects that are important: the Lorentz force v x B on the motion of electrons 
and the coupling of the electron spin to the magnetic field. We first discuss the spin 
coupling and postpone the Lorentz force for the next section. 

An electron can have two spin states, "up" and "down". We'll denote the spin state 
by a discrete label: s = 1 for spin up; s = — 1 for spin down. In a background magnetic 
field B, the kinetic energy picks up an extra term 



where Hb = \e\h/2mc is the Bohr magneton. (It has nothing to do with the chemical 
potential. Do not confuse them!) 

Since spin up and spin down electrons have different energies, their occupation num- 
bers differ. Equation (3.50) is now replaced by 



-E'spin — HbBs 



(3.58) 
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Our interest is in computing the magnetization, which is a measure of how the energy 
of the system responds to a magnetic field, 



M = - (3.59) 

Since each spin has energy (3.58), the magnetization is simply the difference in the in 
the number of up spins and down spins, 

M = /, B (7V T - Nl ) = ^- [h l2 {ze^ B ) - f s/2 (ze-f*» B )] 

At high temperatures, and suitably low magnetic fields, we can approximate / 3 / 2 (z) ~ z 
as z — > 0, so that 

M « ^bVz smh((3fi B B) 

We can eliminate the factor of z in favour of the particle number, 

2Vz 

N = N^ + N^ — cosh((3/2 B B) 

From which we can determine the high temperature magnetization, 

M « /ji B Ntimh((3iji B B) 

This is the same result that you computed in the first problem sheet using the simple 
two-state model. We see that, once again, at high temperatures the quantum statistics 
are irrelevant. The magnetic susceptibility is a measure of how easy it is to magnetise 
a substance. It is defined as 

dM 

X - W (3.60) 

Evaluating x m the limit of vanishing magnetic field, we find that the magnetization 
is inversely proportional to the temperature, 

, (B = o) = M 
AV ; k B T 

This 1/T behaviour is known as Curie's law. (Pierre, not Marie). 
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The above results hold in the high temperature limit. In contrast, at low tempera- 
tures the effects of the Fermi surface become important. We need only keep the leading 
order term in (3.54) to find 

where we're written the final result in the limit hbB <C Ep and expanded to linear order 
in B. We could again replace some of the factors for the electron number N. However, 
we can emphasise the relevant physics slightly better if we write the result above in 
terms of the density of states (3.19). Including the factor of 2 for spin degeneracy, we 
have 

M « fi 2 B g(Ep)B (3.61) 

We see that at low temperatures, the magnetic susceptibility no longer obeys Curie's 
law, but saturates to a constant 

What's happening here is that most electrons are embedded deep within the Fermi 
surface and the Pauli exclusion principle forbids them from flipping their spins in re- 
sponse to a magnetic field. Only those electrons that sit on the Fermi surface itself - 
all g{Ep) of them — have the luxury to save energy by aligning their spins with the 
magnetic field. 

Notice that x > 0- Such materials are called paramagnetic. The effect described 
above is known as Pauli paramagnetism. 

3.6.7 Landau Diamagnetism 

For charged fermions, such as electrons, there is a second consequence of placing them 
in a magnetic field: they move in circles. The Hamiltonian describing a particle of 

— * — * 

charge — e in a magnetic field B = V x A is, 




This has an immediate consequence: a magnetic field cannot affect the statistical 
physics of a classical gas! This follows simply from substituting the above Hamiltonian 
into the classical partition function (2.1). A simple shift of variables, ff = p + eA/c 
shows that the classical partition function does not depend on B. This result (which was 
also mentioned in a slightly different way in the Classical Dynamics course) is known 
as the Bohr-van Leeuwen theorem: it states that there can be no classical magnetism. 
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Let's see how quantum effects change this conclusion. Our first goal is to understand 
the one-particle states for an electron in a constant background magnetic field. This is 
a problem that you will also solve in the Applications of Quantum Mechanics course. 
It is the problem of Landau levels. 

Landau Levels 

We will consider an electron moving in a constant magnetic field pointing in the z 
direction: B = (0, 0, B). There are many different ways to pick A such that B = Vxi 
gives the right magnetic field. These are known as gauge choices. Each choice will give 
rise to the same physics but, usually, very different equations in the intervening steps. 
Here we make the following gauge choice, 



A = (-By, 0,0) 



(3.62) 



Let's place our particle in a cubic box with sides of length L. Solutions to the 
Schrodinger equation Hip = Etp are simply plane waves in the x and z directions, 

But the wavefunction f{y) must satisfy, 

h 2 d 2 



f(y) = E'f(y) 



2m dy 2 

But this is a familiar equation! If we define the cyclotron frequency, 



(3.63) 



eB 
mc 



Then we can rewrite (3.63) as 

h 2 d 2 



f(y) = E'f(y) 



where yo = hk x c/eB. This is simply the Schrodinger equation for a harmonic oscillator, 
jiggling with frequency u c and situated at position y . 

One might wonder why the harmonic oscillator appeared. The cyclotron frequency 
uj c is the frequency at which electrons undergo classical Larmor orbits, so that kind of 
makes sense. But why have these become harmonic oscillators sitting at positions yo 
which depend on k x l The answer is that there is no deep meaning to this at all! it is 

— * 

merely an artifact of our gauge choice (3.62). Had we picked a different A giving rise to 
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the same B, we would find ourselves trying to answer a different question! The lesson 
here is that when dealing with gauge potentials, the intermediate steps of a calculation 
do not always have physical meaning. Only the end result is important. Let's see what 
this is. The energy of the particle is 



E = E' + 



h 2 k 2 z 

2m 



where E 1 is the energy of the harmonic oscillator, 



E' 



n 



huj. 



These discrete energy levels are known as Landau levels. They are highly degenerate. 
To see this, note that k x is quantized in unit of Ak x = 2-n/L. This means that we 
can have a harmonic oscillator located every Ay = 2-Khc/eBL. The number of such 
oscillators that we can fit into the box of size L is L/Ay . This is the degeneracy of 
each level, 



eBL 2 
2nhc 



$ 



where $ = L 2 B is the total flux through the system and $ = 2irhc/e is known as the 
flux quantum. (Note that the result above does not include a factor of 2 for the spin 
degeneracy of the electron). 

Back to the Diamagnetic Story 

We can now compute the grand partition function for non-interacting electrons in a 
magnetic field. Including the factor of 2 from electron spin, we have 



, ^ L f ~ 2L 2 R 
]ogZ = -J dk,^ — ]og 

J n=0 u 



z exp 



ph 2 k 2 

2m 



Pfiw c (n + 1/2) 



To perform the sum, we use the Euler summation formula which states that for any 
function h(x), 

00 roc 

]T h(n + 1/2) = / h(x)dx + -h'(0) + ... 



We'll apply the Euler summation formula to the function, 

(3h 2 k 2 



h(x) = J 



dk z log 



1 + exp 



2m 



+ (3x 
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So the grand partition function becomes 



oo 



log Z = ^ J2 h ^ - ^ n + V2)) 



VBhw c dh{^) 
nu c x) ax — — : h . . . 



Vm 
2tt 2 ^ 2 



. . (3.64) 



e /3(h?k 2 /2m-fi) _|_ J 



The first term above does not depend on £>. In fact it is simply a rather perverse way 
of writing the partition function of a Fermi gas when _B = 0. However, our interest 
here is in the magnetization which is again defined as (3.59). In the grand canonical 
ensemble, this is simply 

M= l51ogZ 

13 8B 

The second term in (3.64) is proportional to B 2 (which is hiding in the ou 2 term). 
Higher order terms in the Euler summation formula will carry higher powers of B and, 
for small B, the expression above will suffice. 

At T = 0, the integrand is simply 1 for \k\ < kp and zero otherwise. To compare 
to Pauli paramagnetism, we express the final result in terms of the Bohr magneton 
Hb = \e\h/2mc. We find 

M = -&g{E F )B 

This is comparable to the Pauli contribution (3.61). But it has opposite sign. Sub- 
stances whose magnetization is opposite to the magnetic field are called diamagnetic. 
The effect that we have derived above is known as Landau diamagnetism. 
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